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摘 要 : 


中 小 型 河流 的 径流 量 精准 监测 对 干旱 区 生态 稳定 具有 重要 意义 。 然 而 中 小 型 河流 流量 还 


TR 


感 精准 反 演 存在 困难 。 以 新 疆 尼 勒 克 县 境内 的 喀什 河 种 峰 场 河 段 为 例 , 基 于 关系 拟 合法 ,依据 实 
测 水文 数 据 、 无 人 机 数据 和 卫星 数据 ,构建 河 宽 水 深 与 流量 之 间 究 函数 关系 模型 ,并 利用 卫星 数据 
的 时 序 性 , 反 演 监 测 河 段 24 次 不 同时 期 的 径流 量 。 反 演 结果 表明 : 当 径流 量 为 0~50 ms! Ft 50~ 
100 ms!+ 时 ,基于 河 宽 的 水 力 儿 何 形态 径流 量 反 演 效 果 最 优 , 均 方 根 误差 (RMSE) 分 别 为 7.15 m s" 


和 2.81 ms.s-; 当 径流 量 为 100~200 ms 和 >200 ms.s 时 ,基于 水 深 和 河 宽 的 水 力 几 何 形态 径流 量 


反 演 效果 最 佳 ,RMSE 分 别 为 13.37 m-s! 1.06 ms.s-!。 研 究 结果 可 为 水 文 资 料 缺 乏 区 的 中 小 型 河 
流 径流 精细 化 监测 与 管理 提供 一 种 新 方法 ,也 对 洪流 灾害 预测 、 水 能 资源 开发 与 水 生态 系统 修复 


具有 较 高 的 参考 价值 。 


I S 词 : 无 人 机 遥感 ; Sentinel-2; 河道 流量 ; 关系 拟 合法 ; 反 演 
文章 编号 : 1000 - 6060(2023)03 - 0385 - 12(0385 ~ 0396) 


径流 量 是 水 资源 开发 利用 河流 生态 修复 与 保 
护 所 需 的 基础 数据 ,也 是 评估 水 资源 储量 能源 风 
险 ,社会 经 济 发 展 的 重要 依据 "”。 但 很 多 地 区 受到 
地 域 环 境 经济 .政策 等 因素 的 限制 ,固定 水 文 观测 
站 点 分 布 稀 玖 其 至 缺失 ,造成 水 文 要 素 信息 缺乏 。 
为 了 解决 在 水 文 无 资料 或 少 资料 地 区 获取 有 效 的 
流量 数据 的 科学 问题 ,很 多 学 者 将 遥感 技术 应 用 于 
河道 的 流量 监测 与 反 演 ,在 基础 理论 与 实际 应 用 方 
面 进行 大 量 研 究 "。 

随 着 地 理 空间 信息 科学 的 发 展 ,基于 关系 拟 合 
法 ,通过 水 力 特 征 变量 反 演 径 流量 的 方法 已 成 为 遥 
感 领域 的 一 个 重要 研究 方向 。 通 过 关系 拟 合 法 来 
耦合 遥感 数据 反 演 径流 量 的 科学 问题 ,国内 外 学 者 
已 进行 大 量 研究 。Leopold 等 "通过 研究 美国 西部 
河流 的 水 文 数据 ,首次 提出 将 河 宽 作 为 水 力 特征 变 
ti ,构建 河 宽 与 流量 的 震 函 数 模型 进行 径流 量 反 
演 ,该 理论 为 径流 量 反 演 提供 了 新 的 思路 ;随后 
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Rhodes ^ 在 全 球 范围 内 选取 587 组 河流 ,以 河 宽 作 
为 水 力 特征 变量 反 演 流量 ,验证 了 关系 拟 合法 的 可 
靠 性 ;在 前 者 的 研究 基础 上 ,Gleason 等 "又 利用 遥 
感 卫 星 数据 获取 长 江 、 密 西西 比 河 、 亚 达 斯 巴 卡 河 
等 大 型 河流 的 河 宽 ,构建 了 多 站 水 力 儿 何 形态 的 函 
数 关系 进行 径流 量 反 演 ,进一步 推广 了 该 方法 的 应 
用 ;Bjerklie 等 ”通过 河 宽 水深 与 流速 等 水 力 特征 
参数 构建 径流 量 反 演 模型 模拟 流量 ;Sichangi 4 
利用 多 源 遥 感 数据 获取 河 宽 和 水 深 等 水 力 特征 变 
量 , 对 尼罗河 和 亚马逊 的 径流 量 进行 反 演 。 基 于 国 
外 的 研究 成 果 , 李 甲 振 等 " 概 化 河道 断面 ,构建 概 
化 断面 的 河 宽 与 流量 的 过 函数 模型 ,通过 实测 河 宽 
对 径流 量 进行 反 演 ; 赵 长 森 等 ”将 无 人 机 航 测 数 据 
应 用 到 径流 量 反 演 中 ,对 济南 小 清河 流域 的 65 组 河 
流 断 面 进行 径流 量 反 演 ; 杨 胜 天 等 通过 无 人 机 航 
空 摄影 获取 高 精度 的 遥感 影像 ,结合 曼 宁 公式 对 新 
疆 典 型 无 资料 区 中 小 型 河流 径流 量 进行 反 演 。 
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虽然 最 新 研究 已 将 低空 冰 感 数据 与 卫星 数据 
相 结 合 ,对 水 文 资料 缺乏 区 的 中 小 型 河流 径流 量 进 
行 反 演 ,但 相关 研究 仍然 需要 进一步 完善 ””。 本 
文 利 用 无 人 机 数据 的 高 精度 优势 与 卫星 数据 的 时 
序 性 特征 ,将 低空 冰 感 .卫星 遥感 和 实测 水 文 数据 
三 者 相 结 合 进行 径流 量 反 演 人 研究 ,利用 河 宽 、 水 深 
和 流量 数据 构建 舌 函 数 关系 模型 ,而 合 多 光谱 遥感 
数据 ,实现 对 中 小 型 河道 的 时 序 性 流量 反 演 。 为 实 
现 水 文 信息 缺乏 区 河道 流量 时 序 性 监测 提供 了 一 
种 高 效 .便捷 的 新 途径 ,也 可 用 于 洪流 灾害 预测 与 
水 资源 精细 化 管理 。 


1 研究 区 概况 


本 文选 取 喀 什 河 上 游 出 山口 处 的 种 蜂 场 水 文 
监测 断面 为 研究 对 象 (图 1)。 喀 什 河 地 处 尼 勒 克 县 
境内 , 北 东 南 三 面 环 山 ,发 源 于 依 连 哈 比尔 条 山 , 流 
癌 从 东 向 西 。 河 流 总 长 度 315 km, 流域 面积 9541 
km? ,支流 多 在 北岸 ,流域 呈 狭 长 的 柳 叶 形 ,流域 高 
程 800~4600 m。 河 流 多 年 平均 流速 为 127.9 mes”, 
多 年 平均 径流 量 为 3.902x10" meso FAK YALA UK 
雪 融 水 为 主 ,降水 补给 次 之 ”。 


>z 


at K'E Cove 


图 例 
A 种 峰 场 控制 断面 
一 喀什 河 


高 程 /m 0 23km 
mx. 4573 — 
EM 697 


图 1 尼 勒 克 段 喀什 河水 系 分 布 
Fig. 1 Distribution of water system of Kashi River in 


Nilka section 
2 数据 与 方法 


2.1 研究 方法 

本 文 基于 关系 拟 合法 ,将 低空 遥感 数据 .实测 
水 文 数据 和 卫星 数据 相 结 合 ,构建 军 函 数 模型 ,对 
河流 流量 进行 反 演 ,具体 步骤 如 图 2 所 示 。 通 过 无 
人 机 航拍 获得 的 数字 表面 模型 (DSM ) 数 据 ,对 种 蜂 


场 河 段 均匀 分 割 , 概 化 出 该 河 段 的 水 力 特征 关系 ; 
基于 河 段 水 力 特征 关系 ,利用 Sentinel-2 遥 感 数据 
确定 流量 反 演 所 需 的 经 验 参数 值 ,进行 河流 流量 反 
演 , 并 对 反 演 结果 进行 精度 验证 2” 。 
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iE: DOM 为 正 射影 像 ;DSM 为 数字 表面 模型 。 
2 水 力 几 何 形态 流量 反 演 流程 图 


Fig. 2 Hydraulic geometry flow over inversion flow chart 


211 水 力 将 征 关系 构建 

(1) 河 段 划分 。 河 道 在 自然 发 育 过 程 中 由 于 水 
流 的 侵蚀 .搬运 .堆积 作用 ,难以 使 用 统一 的 关系 表 
述 不 同位 置 的 断面 形状 特征 。 目 标 河 段 的 划分 产 
生 多 个 较 短 的 河 段 ,在 这 些 较 短 的 河 段 内 可 以 归纳 
出 表征 断面 形状 关系 的 河 宽 -水 深 关系 曲线 。 

(2) 划分 河 段 断面 的 河 宽 -水 深 关 系 曲 线 。 提 
取 划 分 河 段 的 高 程 数 据 , 获 取 不 同 水 深 对 应 的 河 宽 
数据 ,绘制 划分 河 段 断面 的 河 宽 - 水 深 关系 曲线 。 

(3) 拟 合 河 段 水 力 特 征 关系 平均 曲线 方程 。 基 
于 划分 河 段 断面 的 河 宽 -水 深 关 系 曲线 ,计算 相同 
水 深 各 个 分 割 河 段 河 宽 的 平均 值 ,使 用 Origin 软件 
拟 合 表示 目标 河 段 的 断面 形状 特征 的 河 宽 -水 深 关 
系 曲线 ,获取 函数 方程 ,实现 输入 任意 河 宽 ( 水 深 ) 
获取 对 应 水 深 ( 河 宽 ) 的 结果 。 

(4) 流量 反 演 。 基 于 关系 拟 合法 ,构建 流量 与 
河 宽 .水深 的 函数 关系 ,通过 卫星 数据 提取 河 段 的 
水 域 面积 ,利用 式 (4) 获 取 概 化 河 宽 ,进而 计算 相应 
水 深 ;由 于 实测 水 文 数据 没有 记录 实际 河 宽 ,而 本 
文 流量 反 演 方法 中 的 经 验 参 数 (4、B、C) 需 要 河 宽 数 
据 率 定 , 故 将 水 深 代 入 拟 合 方程 , 即 可 获取 河 宽 , 进 
而 率 定 经 验 参 数 。 最 终 通过 观测 河 宽 或 水 深 的 变 
化 即 可 对 河流 流量 进行 反 演 研究 。 

2.1.2 深度 学 习 卷 积 神 经 网 络 模型 ”深度 学 习 卷 积 
神经 网 络 是 一 种 基于 数据 表征 的 机 器 学 习 算 法 。 
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它 利 用 卫星 影像 的 空间 信息 和 光谱 信息 ,深层 次 的 
提取 图 像 特 征 ,而 且 对 遥感 影像 的 水 体 识别 具有 和 较 
高 的 准确 率 , 优 于 BP 神 经 网 络 . 文 持 向 量 机 、 归 一 化 
水 体 指数 (NDWT) 等 方法 ,避免 了 传统 方法 中 复 
杂 的 预 处 理 过 程 。 故 本 文选 取 Keras 卷 积 神 经 网 络 
对 遥感 影像 中 的 水 体 进行 提取 ,获取 水 域 面 积 。 图 
3 为 基于 Keras 卷 积 神经 网 络 流程 示意 网 ,显示 了 水 
体 提 取 的 整个 过 程 。 

本 文 使 用 的 基于 Keras 卷 积 神经 网 络 的 遥感 图 
像 水 体 识 别 方法 ,主要 有 5 个 步骤 :(1) 图 像 预 处 
理 。 利 用 现 有 的 模型 与 方法 ,抑制 或 消除 图 像 成 像 
时 的 各 种 误差 ,(2) 数据 标记 。 选 取 Sentinel-2 中 
10 m 分 辩 率 的 波段 2(Blue) JE Ez 3 (Green ) 和 波段 4 
(Red) 进 行 假 彩色 合成 ,对 合成 图 片 中 的 水 体 进行 
标记 ,并 保存 标记 水 体 的 多 边 形 顶点 ,创建 一 张 与 
原 图 相同 大 小 的 图 像 ,设置 为 白色 , 读 取 多 边 形 内 
的 像素 点 设置 为 黑色 ;标记 后 的 图 像 ,白色 部 分 表 
示 陆 地 ,黑色 部 分 表示 水 体 ,(3) 图 像 切 分 与 建立 训 
练 集 。 为 实现 模型 对 水 体 识 别 ,将 遥感 影像 切 分 为 
16x16 的 小 块 ,结合 标记 的 图 像 建立 训练 集 ,(4) FE 
于 Keras 卷 积 神经 网 络 构建 模型 及 训练 。 本 文 的 
Keras 卷 积 神经 网 络 共 设置 3 个 卷 积 层 (Conv) ,3 个 
最 大 池 化 层 (Max Pool) , 1 个 输入 层 (Input) ,2 个 连 


16x16xK| 


Conv 


256x256 pixels F 


注 : 开 为 通道 数 ;pixels 为 像素 ;Conv HABA 


128x128 pixels 64x64 pixels 64x64 pixels 32x32 pixels 
Max Pool 
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接 层 (F) 和 1 个 输出 层 (Output) 。 其 中 , 卷 积 层 和 最 
大 池 化 层 主要 用 于 遥感 图 像 的 水 体 提 取 。(5) 使 用 
训练 后 的 模型 自动 化 识别 水 体 。 通 过 前 向 传播 和 
反 向 传播 2 个 重要 训练 过 程 ” ,确定 了 该 模型 的 参 
数 , 最 后 将 假 彩 色 遥 感 图 像 输 入 , 即 可 获取 水 体 
言 息 。 

2.1.3 水 力 参 数 计 算 在 众多 水 文学 预测 模型 中 ， 
水 力 梯度 和 焰 率 是 不 可 或 缺 的 重要 参数 ,但 已 有 人 研 
究 证 明 ,在 水 力 梯 度 、 糙 率 缺 失 的 情况 下 依然 可 以 
实现 对 流量 的 精准 预测 ”。 关 系 拟 合法 基于 河 宽 、 
JU 5 TALE AK TI PRAES 3 TAE E PR RS 
现 河 道 径流 量 反 演 。 关 系 拟 合法 是 适用 短 消 数 型 
断面 、 三 角形 断面 河流 和 冤 深 比较 大 的 河流 进行 流 
量 反 演 *。 计 算 公式 如 下 : 


QzAW (1) 
Q- BI (2) 
Q-CWiP (3) 


TA LALBLCOS ZEE C SQ LH WAA AI o 
(m? s) Km) 河 宽 m)。 

由 于 利用 遥感 数据 直接 获取 河 宽 , 可 能 会 受到 
山体 阴影 云层、 冰雪 等 因素 的 干扰 , 河 宽 会 产生 较 
大 的 误差 。 因 此 ,利用 河 段 的 水 域 面积 与 长 度 的 比 


Conv Max Pool 


16x16 pixels 32x32 pixels 


Max Pool Conv 


E: Max Pool Jj WWE s F NEBR. 


图 3 Keras 卷 积 神经 网 络 流 程 图 


Fig. 3 Flow chart of Keras convolutional neural network 
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值 ,产生 的 误差 较 小 ,可 获得 精确 度 较 高 的 的 河 宽 
数据 。 计 算 公式 如 下 


W=S/L (4) 
AP: WS LAP al Ariel GE (m) KERE Om) 、 河 段 


长 度 (m)。 

2.1.4 精度 评估 本 文采 用 3 种 精度 评估 指标 来 验 
证 径流 量 反 演 结果 ,分 别 为 相对 精度 (RA)、 均 方 根 
误差 (RMSE) 和 平均 百分比 误差 (MPE)。 计 算 公式 
如 下 : 


Q. - Q.] 
RA- 5 
0, ws 
RMSE- ze (6) 
1x 00. 
MPE- ly “17 s 7 
E Q. (7) 


SUF: QA dii tat c TRES PRU EE (m * s ;0, 为 实测 径 
流量 值 (ms );n 为 反 演 径 流量 值 的 次 数 。 

2.2 数据 来 源 与 处 理 

2.2.1 X. £ i£ Eh d& Sentinel-2 遥 感 数 据 的 空间 
分 辩 率 较 高 ,周期 短 , 且 预 处 理 方 法 比较 简易 ,特别 
是 在 水 域 识别 与 分 析 方 面 被 普遍 应 用 。 本 研究 采 
JH T 2016—2019 年 24 次 不 同时 期 的 Sentinel-2 遥 
感 影像 ,与 实测 数据 的 时 间 相对 应 ,基于 GEE 平 台 ， 
根据 Sentinel-2 If] QA60 波段 特征 ,设置 冰 值 ,筛选 
人 研究 区 域内 的 遥感 数据 。 选 用 Sentinel-2 中 10 m 分 
辩 率 的 波段 2(Blue) 波段 3(Green) 和 波段 4(Red ) , 
进行 假 彩 色 合 成 ,使 用 Keras 卷 积 神经 网 络 模型 进 
行 水 体 提取 。 

2.2.2 低 室 遥感 数据 ”本文 无 人 机 数据 取 自 于 2018 
年 4 月 尼 勒 克 喀 什 河 种 蜂 场 河 段 ,该 时 间 段 喀什 河 处 
于 枯水期 ,水 位 低 流速 组 ,水深 一 般 不 超过 50 cm, 
可 以 最 大 程度 上 获取 河道 地 形 数据 。 本 研究 使 用 
Keg Phantom 4 Advanced 获 取 高 精度 的 低空 无 人 机 
地 形 数据 ,选择 DJI GS Pro 软件 控制 无 人 机 飞行 , 规 
划 飞 行路 线 与 航 测 任务 。 获 取 种 蜂 场 河 段 的 高 精 
度 地 形 数据 主要 有 以 下 几 步 :(1) 设置 无 人 机 飞行 
轨迹 。 依 据 种 蜂 场 水 文 监测 断面 的 地 理 位 置 ,飞行 
高 度 一 般 为 100 m, 以 河道 为 中 心 设置 飞行 航线 。 
为 了 保证 多 次 飞行 时 起 飞 点 的 一 致 性 ,应 用 点 标记 
无 人 机 落 点 的 位 置 坐 标 及 其 周边 参考 地 物 的 位 
置 。(2) 无 人 机 空中 作业 。 无 人 机 每 次 飞行 的 高 度 、 
航线 要 确保 一 致 ,各 航 带 距离 地 面 高 度 也 要 一 致 , 确 


保 邻 近 各 景 无 人 机 影像 的 重 羡 率 为 70%~80% ,方便 
后 续 影 像 的 几何 校正 和 灸 能 。(3) 无 人 机 影像 拼接 
与 处 理 。 利 用 Pix4D mapper 软 件 自动 化 处 理 无 人 机 
影像 ,生成 喀什 河 种 蜂 场 河 段 的 正 射 影像 (DOM)、 
数字 表面 模型 (DSM)。 最 后 在 ArcGIS 中 的 3DAna- 
lyst 模 块 上 ,对 种 蜂 场 河 段 所 在 的 200 m 范 围 内 DSM 
影像 均匀 提取 4 个 断面 地 形 数据 ,再 将 水 上 断面 地 
数 据 与 测 杆 实测 水 深 数据 相 结 合 , 生 成 完整 的 河 
道 断面 。 

2.2.3 实测 数据 本 文 实测 数据 来 源 于 种 蜂 场 水 文 
站 2009 一 2019 年 的 实测 流量 和 水 位 数据 (图 4)。 为 
了 验证 流量 反 演 的 可 靠 性 ,基于 河 宽 与 水 深 之 间 的 
水 力 几 何 模型 ,将 实测 水 文 数据 分 为 两 部 分 , 率 定 
期 (2009 一 2016 年 ) 与 验证 期 (2016 一 2019 年 ) , 率 定 
期 的 数据 被 用 于 率 定 式 (1)~ 式 (3) 中 的 4.B、C3 个 
经 验 参数 ,验证 期 则 用 于 反 演 径流 量 ,并 对 反 演 结 
果 进 行 精度 评价 。 
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图 4 种 蜂 场 水 文 站 径流 量 -水 位 关系 曲线 


Fig.4 Runoff-water level relationship curve of 


Zhongfengchang hydrological station 


3 结果 与 分 析 


3.1 河道 断面 与 水 深 关系 曲线 

种 蜂 场 水 文 站 所 在 的 200 m 河 段 均匀 的 提取 4 
个 断面 ,利用 ArcGIS 中 的 3DAnalyst 模 块 在 DSM 影 
像 获取 分 割断 面 的 高 程 ,绘制 每 个 断面 的 河 宽 -水 
深 关 系 曲 线 (图 5)。 
3.2 拟 合 河 段 水 力 特征 关系 平均 曲线 

基于 种 蜂 场 河 段 4 个 断面 的 水 力 特 征 关 系 曲 
线 , 计 算 该 河 段 4 个 断面 同一 水 深 的 平均 河 宽 ,绘制 
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图 5 河道 断面 与 水 力 特征 关系 曲线 


Fig. 5 Relationship curve between river cross-section and hydraulic characteristics 


河 宽 - 水 深水 力 特征 关系 的 平均 曲线 ,作为 目标 河 
段 的 断面 形态 特征 (图 6)。 

本 文 使 用 Origin 软件 获取 河 宽 - 水 深水 力 特 征 
关系 曲线 的 拟 合 方程 。 该 数值 分 析 软 件 根 据 麦 考 
特 迭 代 法 (Levenberg-Marquardt) 的 全 局 优化 算法 对 
数据 进行 迭代 计算 ,计算 出 效果 最 佳 的 拟 合 方程 及 
残 差 分 析 ( 表 1)。 由 表 1 可知 ,图 6a 与 图 6b 的 曲线 
拟 合 殖 果 较 优 ,相关 系数 均 达 到 98% , 残 差 平方 和 
分 别 为 4.77、3.51, 均 方 根 误差 分 别 为 0.56、0.50, 准 
确 反 映 了 种 蜂 场 河 段 河 宽 与 水 深 之 间 的 变化 特 
征 。 根 据 表 1 中 下 fF( 风 方程 ,利用 式 (4) 获 取 概 化 
河 宽 , 即 可 计算 对 应 水 深 。 由 于 实测 水 文 数据 中 仅 
记录 了 水 深水 位 和 径流 量 数 据 , 而 没有 记录 河 宽 
数据 ,但 可 以 通过 =C( 四 方程 可 获取 河 宽 数 据 , 完 
成 关系 拟 合 法 中 率 定期 需要 河 宽 参 与 经 验 参 数 4、 
B.C 的 率 定 。 

3.3 水 力 特 征 参数 提取 

基于 Sentinel-2 卫星 数据 ,通过 深度 学 习 Keras 
卷 积 神经 网 络 模 型 获取 长 度 为 200 m 的 种 蜂 场 河 段 
的 水 域 面积 。 在 获取 水 域 面积 时 , 因 无 人 机 遥感 数 
据 与 Sentinel-2 卫星 遥感 数据 并 未 重合 , 故 以 Senti- 
nel-2 卫星 遥感 数据 为 基准 ,上 .下 游 边 界 各 扩展 1 
个 像素 宽 (10 m) ,进而 获取 目标 河 段 的 概 化 河 宽 。 
该 Keras 卷 积 神经 网 络 模型 提取 准确 率 高 达 
99.54% , 误 报 率 为 0.76%。 表 2 显示 了 种 蜂 场 河 段 
24 次 不 同时 期 的 水 域 面 积 、 概 化 河 宽 和 概 化 水 深 , 通 
过 式 (4) 计 算得 到 的 概 化 河 宽 , 利 用 表 1 中 天 所 有 取 方 
程 即 可 获得 概 化 水 深 。 

3.4 流量 反 演 
基于 关系 拟 合法 ,根据 率 定期 实测 水 文 数据 来 


确定 经 验 参 数 4、B、C, 实 现 对 各 级 流量 反 演 。 如 图 
7 所 示 , 图 7a~c 反 映 了 各 公式 对 各 级 流量 的 反 演 效 
HE ,图 74 则 是 对 比 各 公式 对 各 级 流量 的 反 演 效 
果 。 结 合 表 3 可 看 出 , 当 径 流量 0~50 ms fI50- 
100 ms 时 ,以 河 宽 作 为 水 力 特 征 变量 的 式 (1) 反 演 
结果 最 佳 ,RMSE 分 别 为 7.1$ ms 2.81 ms; 44% 
流量 >200m.s 时 ,以 河 宽 , 水深 作为 水 力 特征 变量 
的 式 (3) 反 演 结 果 最 优 ,RMSE 为 1.06 m? s ,与 实测 
数据 基本 吻合 ; 当 径流 量 100-200 m- sE}, LLKR 
为 水 力 特 征 变量 的 式 (2),RMSE X 13.60 m-s”, ifi; 
(3) 的 RMSE 为 13.37 ms.s…, 式 (2) 反 演 结果 略 差 于 
式 (3)。 整 体 可 以 看 出 , 当 径 流量 在 0~100m.s FT, 
与 式 (1) 相 比 , 式 (2) 反 演 结果 整体 误差 较 大 ;而 径 
流量 >100 m? «s "时 ,整体 误差 较 少 ,但 反 演 精度 略 差 
于 式 (3)。 


4 讨论 


4.1 研究 结果 对 比分 析 

近年 来 ,许多 学 者 通过 各 种 遥感 数据 对 径流 量 
进行 反 演 , 并 取得 了 重大 进展 。 王 鹏 飞 等 ”在 水 力 
几何 形态 径流 量 反 演 研 究 中 ,使 用 的 研究 方法 与 本 
文 一 致 ,构建 径 流量 与 河 宽 , 水深 的 需 函 数 模型 ,对 
径流 量 进行 反 演 , 并 进行 精度 评价 ,结果 发 现 :径流 
量 <200 ms 时 , 式 (1) 模 拟 效 果 较 好 ,RMSE 为 
4.65 m? -s~ ; 当 径 流量 >200 ms 时 , 式 (3) 的 精度 评 
价 结果 较 好 ,平均 合格 率 为 80% ,但 其 反 演 精度 仍 
可 进一步 提高 。 在 本 研究 中 ,将 验证 期 实测 流量 数 
值 (最 小 值 23.20 m’…s', 最 大 值 367 m8) ,详细 
分 为 4 个 等 级 :=>200 m +s"! 100~200 m^: s, 50~ 
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250 - (d) 断面 2 (水 深 与 河 宽 关 系 ) 
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图 6 种 蜂 场 河 段 4 个 断面 的 拟 合 水 力 特征 曲线 


Fig.6 Fitting hydraulic characteristic curves of four sections at Zhongfengchang river section 


100 m° +s! 和 0~50m.s 。 由 表 3 可 知 : 当 径流 量 
«100 ms! , 式 (1) 的 各 项 评价 指标 均 优 于 式 (2) 和 式 
(3) ,其 平均 合格 率 为 90%。RMSE .MPE 2 个 评价 指 
标 均 表明 径流 量 50~100 ms 时 , 反 演 结果 最 佳 , 反 


演 精 度 优 于 同类 研究 的 结果 ; 当 径 流量 0~50 ms" 
时 , 反 演 效果 次 之 ,对 水 资源 精细 化 管理 有 一 定 
的 参考 价值 ; 当 径 流量 100~200 n? s” AY, È (3) 
演 效果 较 好 ,整体 优 于 其 他 公式 ,RMSE 均值 为 
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于 山区 和 平原 的 河流 断面 


BAS (Al 8e) ,表现 为 一 边 


Tab. 1 Morphological characteristic equations and fitting 


curves of the generalized river section 


函数 H=F(W) W=G(H) 

方程 y=BotB + Box? + By +Bax* y=BotB xt Bux Ba Bax? 

参数 值 Bo:37.5882;B::30.8797 —— B,:88.38032; B, :3.59985 
B.:-3.65321;B.:0.36687 — B2:—0.05631;B3:0.00429 

B.:-0.0143 B.:-0.000002 

相关 系数 0.98 0.98 

均 方 根 误差 0.56 0.50 

残 差 平方 和 4.77 3.51 


TE: FW ARE WAY PRL GN AEE H A RRG HAKR SW 2g 
His Bo Bi B2 Bs BABE y Tg 因 变 量 ;x 为 自 变 量 。 


13.37 n? s , MPE 为 -0.05, 反 演 结果 的 误差 均 小 于 
实际 流量 的 20%; 当 径 流量 >200 m-s 1! 时 , 式 (3) 反 
演 结果 最 优 ,RMSE 为 1.06 m-s, RA 290.3896, M 
体 而 言 , 当 种 蜂 场 河 段 径流 量 较 高 时 ,把 水 深 和 河 
宽 作 为 输入 变量 , 式 (3) 反 演 结果 整体 较 优 。 
4.2 典型 河道 的 适用 性 

喀什 河 种 蜂 场 水 文 断面 是 典型 的 宽 深 比 较 大 
的 三 角形 断面 ,其 河 宽 水 深 随 径流 波动 而 变化 ,可 
较 好 地 表征 径流 量变 化 。 如 图 8 所 示 的 5 种 典型 断 
面 ,以 河 宽 表征 径流 量变 化 时 ,本 文 的 研究 方法 具 
有 一 定 的 参考 价值 。 三 角形 河道 断面 具有 明显 的 V 
JE Ta FE ERE ,为 典型 山区 河流 断面 形态 ,可 通过 河 
宽 的 波动 表征 径流 的 变化 过 程 (图 8a) ; 宽 浅 型 辫 状 
河道 断面 具有 明显 的 突起 , 河 宽 同样 可 以 表征 径流 
变化 ”( 图 8b); 特 殊 三 角形 河流 断面 是 普遍 存在 


河岸 坡度 较 陡 , 另 一 边 河岸 坡度 较 缓 , 当 河 宽 变 化 
时 ,流量 也 表现 出 与 河 宽 相 同 的 变化 趋势 ;在 平原 
地 区 与 丘陵 地 区 ,和 抢 形 河道 断面 和 复合 型 河道 断面 
为 常见 的 河流 断面 形态 ,表现 为 水 力 梯度 小 ,流速 
缓 (图 8d~e)。 其 中 ,复合 型 河道 断面 宽 深 比较 大 ， 
适用 于 河 宽 表 征 径流 的 变化 ,而 矩形 河道 断面 宽 深 
比较 小 ,用 河 宽 表征 径流 变化 不 适用 于 此 类 断面 。 
本 文 构建 河 宽 -水 深水 力 特征 关系 曲线 ,利用 Senti- 
nel-2 卫星 数据 时 序 性 ,对 24 次 不 同时 期 的 径流 量 
进行 反 演 , 实 现 径流 量 在 时 空 上 的 扩展 。 该 研究 方 
法 对 河 宽 可 表征 径流 变化 的 河流 ,具有 一 定 的 适用 
性 ;但 该 方法 也 具有 一 定 的 局 限 性 ,不 适用 于 宽 深 
比较 小 的 梯形 断面 河流 或 矩形 断面 的 河流 。 
4.3 影响 因素 及 敏感 性 分 析 

影响 径流 量 反 演 精度 的 因素 主要 包括 河 段 长 
度 与 反 演 过 程 。 首 先是 河 段 长 度 对 流量 反 演 有 着 
较 大 的 影响 , 概 化 河 宽 是 通过 卫星 遥感 影像 提取 水 
域 面积 ,利用 式 (4) 计 算 河 宽 ,但 在 水 域 面积 提取 的 
过 程 中 , 河 段 过 长 必然 使 提取 的 水 域 面积 与 实际 面 
积存 在 一 定 的 误差 ,影响 最 终 的 反 演 结 果 。 所 以 ， 
当 河 段 较 短 时 ,河道 水 文 断面 的 水 力 几 何 特征 更 为 
精确 。 本 文 基于 DSM 数 据 , 在 长 度 为 200 m fr] inf Ez 
中 ,均匀 提取 4 个 断面 ,研究 径流 量 反 演 , 最 终 流 量 反 
演 精 度 验证 了 河 段 长 度 越 短 , 反 演 结果 越 精确 2 。 

其 次 是 流量 反 演 过 程 , 反 演 过 程 中 需要 河 宽 与 
水 深 数据 来 率 定 经 验 参数 ,经 验 参 数 也 会 对 反 演 的 


表 2 喀什 河 种 蜂 场 断面 的 水 域 面积 \ 概 化 宽度 及 水 深 


Tab.2 Watershed area, generalized width and water depth of the Zhongfengchang cross-section in the Kashi River 


日 期 (年 -月 -日 ) 概 化 水 深 /m 概 化 河 宽 /m 面积 /mY 日 期 (年 -月 -日 ) ” 概 化 水 深 /m ” 概 化 河 宽 /m 面积 /7 
2016-09-11 0.91 51,32 10264 2017-10-03 0.81 47.27 9454 
2016-10-09 0.80 47.08 9416 2018-03-08 0.71 43.34 8668 
2016-10-15 0.78 46.06 9212 2018-03-18 0.73 44.03 8806 
2017-03-21 0.70 43.56 8712 2018-04-17 0.79 46.80 9360 
2017-04-10 0.78 46.12 9224 2018-04-25 0.84 48.80 9760 
2017-05-20 0.92 51.99 10398 2018-07-19 1.78 79.83 15966 
2017-06-09 1.70 77.78 15556 2018-08-18 1.83 8124 16248 
2017-07-06 1.90 83.45 16690 2018-08-28 1.74 78.75 15750 
2017-08-08 2.09 88.41 17682 2018-09-04 1.75 78.78 15756 
2017-08-26 1.71 71.89 15578 2018-10-07 0.75 45.30 9060 
2017-09-10 1:55 73.11 14622 2019-03-23 0.73 44.16 8832 
2017-09-17 1.08 57.34 11468 2019-04-05 0.72 43.83 8766 
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图 7 径流 量 反 演 结果 与 实测 数据 
Fig. 7 Runoff inversion results and measured data 
3 分 级 精度 评价 结果 
Tab.3 Evaluation results of grading accuracy 
3k 径流 量 0~50 m-s” 径流 量 50~100 m-s” 

à RA/% RMSE/m°- s” MPE RA/% RMSE/na .sn MPE 
式 (1) 4.35 7.15 -0.15 4.38 2.81 -0.01 
式 (2) 42.34 19.75 0.31 16.32 10.85 -0.01 
式 (3) 19.42 8.86 -0.14 8.89 5.96 -0.20 
径流 量 100~200 m°- s” 径流 量 >200 m-s” 

RA/% RMSE/n* s” MPE RA/% RMSE/m°- s” MPE 
X) 13.60 22.29 -0.05 8.53 28.21 -0.01 
式 (2) 8.71 13.60 -0.10 3.54 11.20 0.00 
式 (3) 8.99 13.37 -0.05 0.38 1.06 -0.01 


注 : RA 为 相对 精度 ; RMSE 为 均 方 根 误差 ; MPE 为 平均 百分比 误差 。 


最 终结 果 产 生 一 定 程 度 上 的 影响 。 如 采 水 文 断面 
率 定期 的 径流 量 数值 变化 幅度 较 大 ,直接 率 定 经 验 
参数 ,可 能 导致 反 演 结果 产生 较 大 的 偏差 。 而 本 文 
采用 流量 分 级 的 方法 ,对 各 级 流量 分 别 率 定 综合 参 
数 ,精确 地 模拟 各 级 流量 ,降低 反 演 结果 的 误差 ,最 
终 模拟 效果 较 好 ,说 明 流量 分 级 的 方法 对 客 函 数 模 
型 反 演 径流 量 有 一 定 的 参考 意义 。 


4.4 断面 地 形 及 水 力 几 何 特征 获取 

内 外 许多 学 者 依据 关系 拟 合 法 ,获取 目标 河 
段 的 水 力 特 征 变量 ,对 径流 量 进 行 反 演 。 目 标 河 段 
的 断面 通常 为 抛物 线 型 FF nemo BARR 
拟 合法 不 适用 于 平原 地 区 宽 深 比较 低 的 梯形 河 段 ， 
但 关系 拟 合法 仍 适用 于 大 部 分 河流 ,特别 是 山区 河 
流 , 或 河流 的 中 、 上 游 河 段 , 这 些 河 段 断面 宽 深 比 通 
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(a) 三 角形 河道 断面 (b) 立 状 河道 断面 


35] LEE: 低空 遥感 结合 ] 
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(c) 特殊 三 角形 河道 断面 


图 8 不 同形 态 的 河道 过 流 断 面 


Fig. 8 Different forms of river overflow cross-sections 


T Be ey, MXE EEK fw m DC. ARC 
的 研究 方法 为 实现 水 文 信息 缺乏 区 河道 流量 时 序 
性 预测 提供 了 一 种 快速 便捷 的 新 途径 ,可 用 于 潜 
流 灾害 预测 与 水 资源 精细 化 管理 。 

在 径流 量 反 演 研 究 中 ,水文 监 测 断 面 的 高 程 数 
据 是 重要 的 实测 水 文 数据 ””。 本 文 在 种 蜂 场 河 段 
干 润 期 ,基于 无 人 机 DSM 数据 ,获取 水 下 地 形 生 成 
整个 河道 的 断面 。 但 是 ,在 河道 常年 不 断 流 或 河道 
水 量 庙 急 的 情况 下 ,获取 其 水 下 断面 地 形 数据 , 进 
行径 流量 反 演 ,可 以 采取 以 下 措施 :(1) 涉 水 测量 获 
取 断 面 地 形 。 当 水 深 较 浅 且 不 超过 50 cm\ 流 速 较 
组 时 ,可 以 通过 测 杆 实测 或 目 视 判 别 水 深 数 据 ,与 
无 人 机 获取 的 水 上 地 形 数据 相 结合 ,生成 完整 的 水 
下 断面 。(2) 声 纳 探测 与 人 工 船 测 获 取 断 面 地 形 。 
当 水 深 较 大 、 河 流 注 急 时 ,可 通过 无 人 船 声 纳 探测 、 
人 工 船 测 结合 测 杆 获 取水 下 断面 数据 。 一 般 情况 
下 声 纳 探 测 适用 于 水 深 超过 50 cm 的 河流 ,而 人 工 
妥 测 量 适 用 于 宽 深 比较 大 的 大 型 河流 ,在 复杂 河流 
断面 的 背景 下 ,可 将 两 者 结合 ,测量 效果 较 优 。 将 
实测 水 深 数据 与 水 面 结 合 即 可 获得 河道 断面 地 
形 。(3) 拟 合 水 面 地 形 数据 获取 河道 断面 。 当 河流 
地 理 位 置 偏远 或 在 山区 时 ,不 易 获 取 实 测 水 下 地 形 
言 息 ,可 以 通过 无 人 机 航 测 获取 的 水 面 地 形 数据 对 
水 下 地 形 拟 合 (图 8) , 拟 合 类 型 可 分 为 抛物 线 .三 角 
形 、 梯 形 等 类 型 ,具体 拟 合 类 型 的 选取 ,可 根据 
河流 的 地 理 环境 与 实际 情况 来 确定 ,(4) 基于 基础 
流量 确定 河流 断面 地 形 。 当 河流 较 宽 ,流量 较 组 
时 ,但 河道 环境 较为 复杂 ,如 立 状 河 ,可 通过 选取 多 
个 不 同 水 次 的 基础 流量 值 (Q) ,结合 实际 情况 ,基于 
某 一 基础 流量 值 ,获取 河道 断面 地 形 ,对 河流 流量 


进行 反 演 (图 9)。 
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注 :0 为 某 一 基础 流量 值 。 
图 9 水 下 断面 拟 合 


Fig.9 Underwater Section fitting 


5 结论 


<i 


本 文 基于 关系 拟 合 法 ,均匀 分 割 河 段 ,构建 目 
标 河 段 的 河 宽 与 水 深 的 函数 方程 , 反 演 径流 量 , 研 
究 河 段 分 割断 面 数量 对 流量 反 演 的 影响 ,分析 流量 
分 级 与 反 演 精度 之 间 的 关系 。 得 到 以 下 主要 结论 : 

(1) 基于 实测 水 文 数据 , 利 用 低空 遥感 数据 和 
卫星 数据 , 拟 合 河 宽 , 水 深 与 流量 的 关系 模型 ,构建 
河 宽 -水 深 平均 关系 曲线 , 概 化 河 段 形状 特征 , 河 段 
分 割 得 越 短 ,其 拟 合 效 果 越 精确 。 震 函数 关系 模型 
也 在 种 蜂 场 河 段 表 现 出 较 好 的 模拟 效果 ,对 反 演 中 
小 型 河流 径流 量具 有 较 高 的 可 靠 性 ,为 无 资料 或 缺 
失 资 料 地 区 的 径流 反 演 提供 了 新 的 解决 方法 。 

(2) 从 流量 分 级 与 反 演 精 度 的 关系 可 看 出 , 当 
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径流 量 >200 m-s 时 ,以 河 宽 水深 作 为 水 力 特 征 变 
量 的 径流 量 反 演 效果 最 优 , 平 均 合 格 率 为 99% ,可 
实现 对 洪流 灾害 的 预测 ; 当 径 流量 100~200 m s" 
时 ,RMSE 为 13.37 m-s” , 反 演 效果 略 差 ; 当 流 量 50~ 
100 ms 时 ,以 河 宽 作 为 水 力 特征 变量 ,径流 量 反 
演 效 果 较 好 ,RMSE X 2.81 m-s”, RA X 4.38% ,利于 
实现 水 资源 精细 化 管理 ; 当 径 流量 <50 mss IB, 
RMSE Jj 7.15 ms ,与 径流 量 50~100 m^- s" FALL, 
反 演 效果 略 差 。 

(3) 本 文 研究 方法 依据 关系 拟 合法 ,利用 高 精 
度 的 低空 遥感 数据 和 长 时 序 的 卫星 遥感 数据 ,实现 
水 文 信 息 缺 乏 区 的 径流 量 的 历史 回溯 ,分 析 研 究 区 
域 径流 量 的 时 空 变化 。 该 方法 也 有 一 定 的 局 限 性 ， 
较 适 用 于 山区 河流 或 宽 深 比 大 的 河流 ,应 用 于 复合 
型 断面 的 河流 ,可 能 会 产生 较 大 的 误差 , 仍 需 进 一 
步 的 研究 。 
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Runoff estimation with low altitude remote sensing and satellite images 


JIANG Leipeng"", DING Jianli^", BAO Qingling"", GE Xiangyu"", 
LIU Jingming'*’, WANG Jinjie™ 
(1. College of Geography and Remote Sensing Sciences, Xinjiang University, Urumqi 830046, Xinjiang, China; 2. Xinjiang Key 
Laboratory of Oasis Ecology, Xinjiang University, Urumqi 830046, Xinjiang, China; 3. Key Laboratory of Smart City and 
Environment Modelling of Higher Education Institute, Urumqi 830046, Xinjiang, China) 


Abstract: Accurate monitoring of runoff from small and medium-sized rivers is of great significance for ecologi- 
cal stability in arid areas. However, it is difficult to accurately retrieve the flow of small and medium-sized rivers 
by remote sensing. Taking the Zhongfengchang river section of Kashi River in Nilka County, Xinjiang, China, as 
an example, this study constructed a power function relationship model between river width, water depth, and dis- 
charge based on the relationship fitting method and measured hydrological data, unmanned aerial vehicle data, 
and satellite data. Using the time series of satellite data, the runoff volume of the monitored river section was in- 
ferred 24 times in different periods. The results show that when the runoff rate is 0—50 m’+s ' and 50-100 m.s, 
the inversion of the runoff rate based on the hydraulic geometry of the river width is optimal, with root mean 
square errors (RMSEs) of 7.15 m’+s ' and 2.81 m’+s ', respectively; when the runoff rate is 100-200 m'*s ' and > 
200 m'*s ', the inversion of the hydraulic geometry based on water depth and river width is the best, with RMSEs 
of 13.37 m'-s ' and 1.06 m°- s`’, respectively. These findings provide a new method for the fine monitoring and 
management of runoff of small and medium-sized rivers in areas lacking hydrologic data and have high reference 
value for flood disaster prediction, hydropower resource development, and water ecosystem restoration. 
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